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TECHNICAL MEMORANDUM X- 64698 


THE PSEUDONOISE TEST SET: COMMUNICATION 
SYSTEM'S PERFORMANCE EVALUATION BASED 
UPON RMS ERROR TESTING 


I. INTRODUCTION 


Until recently, the usual practice in end-to-end testing of communica- 
tion systems has been to feed simple, deterministic signals into the system 
being tested and note the response. The response was then compared with a 
delayed and amplitude scaled version of the input signal to enable one to 
determine the distortion caused by the system and, thus, establish a system 
figure of merit. Unfortunately, most communication systems experience 
complex input signals and true simulations of system inputs are difficult to 
obtain. -Even with the approximation of these complex signals (for instance, 
the sum of a number of independent sine waves) , the deterministic approach 
has often proven unacceptable. Usually the number of component signals 
were inadequate to truly simulate the normal system input. Asa result, 
beat-frequencies or harmonics effects were often missed and the true system 
performance was not evaluated. To overcome this deficiency, band-limited 
noise (usually limited to the same operational bandwidth as the tested system) 
is now ofteji used as a test signal and represents a very general approach to 
communication system performance evaluation. Yet, while band-limited 
noise is an excellent tool with which to simulate realistic data conditions, it 
has proven difficult to uniformly delay — thereby making valid system input/ 
output comparisons difficult to accomplish. 

There are two basic ways to test a system utilizing band-limited noise 
as a test signal — digital and analog (reference Appendix A) . With the 
digital approach, band-limited noise is fed into a test system and the output 
monitored. The system input is sampled and digitized (thereby the digital 
notation) . These digitized samples are compared with digitized samples of 
the test system output yielding mean square values of the distortion introduced. 
Constant gain and average delay factors are compensated so that they do not 
contribute to the distortion introduced by the system being tested since these 
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factors are not usually considered as errors. There are two primary advan- 
tages of such a scheme: A complex test signal is used that is similar to the 
expected operational input spectrum, and all operations after sampling are 
digital. The primary disadvantages are that the effects of system spikes, 
such as FM "clicks,” are often still missed. 

In the analog system approach, band-limited noise is fed into a test 
system. The system output is compared on an rms basis with a delayed and 
scaled version of the system input. The obvious advantages are that the 
equipment is compact and relatively simple to build. The largest disadvantage 
is the difficulty in delaying the complex input without distortion. 

The pseudonoise (PN) test set is designed to provide a reliable, 
operationally simple unit which will allow the previously mentioned analog 
rms end-to-end error measurement of most communication systems to be 
easily performed. It also provides a band-limited pseudorandom noise as 
an input, thereby escaping the disadvantages of most deterministic signals 
while retaining the advantages of using true band-limited noise. In addition, 
the PN test set can be used as a means with which the autocorrelation function 
and impulse response may be determined for certain systems. 

II. THEORY OF PSEUDONOISE TESTING 


Correctly defined noise signals give a good approximation of the data 
signals handled by communication systems. The noise signals must be 
stationary to assure that the measurements are repeatable. Noise can be 
defined by two unrelated parameters, power spectral density and the amplitude 
probability density function. 

Power spectral density (PSD) describes how the energy in a noise 
signal is distributed in frequency. Generally, the noise power spectral 
density, as a function of frequency, can take almost any shape; however, the 
most common specification is that the power spectral density be constant 
with frequency. This describes "white" noise. 

Probability density function (PDF) of a noise signal is defined as that 
function which, when integrated over an interval (Xj, X 2 ) , gives the proba- 
bility that the noise amplitude, X^, lies in the interval (Xj, X^ . Specifically, 
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( 1 ) 


X 2 

P(x t == X n < X 2 ) = / p(x)dx , 


and p(x) is the amplitude probability density function, PDF. 

True noise data signals are difficult to delay without introducing 
appreciable distortion. This problem is eliminated by using pseudorandom 
noise signals. PN signals produced by the test set are, in fact, not random 
but are completely definitive and repetitive. These PN signals can be made 
to closely resemble true random noise signals from both the power spectral 
density and probability density function aspects. 

A pseudorandom binary sequence can be generated by a shift register 
with an input dependent on feedback from certain register stages. The 

maximal length (ML) PN binary sequence is of length L and contains 2° - 1 
bits before repetition of the sequence (where N is the number of register 
stages) . 

An important characteristic of pseudorandom sequences is the corre- 
lation property. If a binary sequence is pseudorandom, the correlation 
function of the sequence, and the sequence shifted with respect to itself, is 
constant for all shifts other than shifts of an integer number of sequence 
periods. Specifically, if vector A is a sequence (a t , a 2 , . . . , a ) , and B 

is a shifted version of this same sequence (with j shifts) , then the corre- 
lation function ( R) is: 


, L 

R (j) = - l 2 < a i © V ( 2 ) 

i= 1 

where the (+) indicates modulo-2 addition. If the sequence is pseudorandom, 
•R(j = 0) = 1 , 

R(0 < j < N) = 1/L , (3) 

R(j = L) = 1 where L = 2^ - 1 
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This correlation property of a pseudorandom sequence is particularly impor- 
tant because it verifies that the signal generation approximates a random 
process. 

The power spectrum envelope of a PN sequence is a (sin X/X) 2 curve 
with the first zero occurring at the sequence clock frequency (co ) . If the 

PN sequence is lowpass filtered, with the cutoff frequency (co^) restrained 

to co < co, /15 , then the PDF of the resultant analog signal can be approxi- 
o k 

mately Gaussian. However, in section VII a situation is described where the 
probability density function of the filtered sequence may not be Gaussian. 


Noise signals generated by lowpass filtering of ML PN sequences are 
partially described by the cutoff frequency and order of the filter. The spec- 
trum amplitude is approximately flat for co < co q , the envelope is down three 

decibels (dB) at co q . At frequencies greater than co q , the envelope rolls off 

at a rate of 6 m dB /octave (where m is the data filter order) . The spectrum 


is not continuous but consists of closely spaced spectral lines 



Figure 1 is a functional block diagram of the PN test arrangement. 
The filtered PN sequence is represented as a sum of sinusoids, 


CO 

Y C cos (cot) . (4) 

i=l 


The filtered sequence is passed through the system with transfer 
function H(s) , and the system output is 

CO 

Y A. cos (co. t + o.) ( 5 ) 

i=l 1 1 1 


A scaled and delayed version of the original analog signal is 


CO 

Y B. cos (co. t + TCO.) . (6) 

i=l 1 
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Figure 1. Block diagram of PN test arrangement. 


The difference between these two complex waveforms is evaluated on an rms 
basis by the PN test set to give a figure of merit relating to the system's 
ability to transfer the complex pseudonoise data signals. 

The PN test set was designed to perform the measurement of the wave- 
form distortion as described above. It has been successfully used to measure 
single-channel and multi-channel performance in CBW/FM communication 
systems. Specific test procedures were outlined to allow for measurement 
of the total end-to-end system error as well as the major individual sources 
of error, namely, channel filter distortion, adjacent channel crosstalk, 
intermodulation, and RF noise. The measured data have been utilized, in 
conjunction with analytically derived data, in the analysis of existing telem- 
etry equipment as well as in the designs proposed for new and improved 
systems. 


III. DESCRIPTION OF A PN TEST SET 


The PN test set as designed and fabricated provides both real-time 
and record/playback systems testing capabilities for measuring rms error, 
system time delay, correlation function, and system impulse response 
(assuming a linear system). 
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Components 


Basically, the PN test set consists of a pseudorandom binary sequence 
generator, sequence delay circuitry, spectrum shaping filters, signal con- 
ditioning amplifiers, a difference amplifier, a system clock, and tape 
recorder interface circuitry. Figure 2 is a simplified block diagram of the 
PN test set. Figure 3 is a photograph of the unit. 

Clock. The system clock, a free-running multivibrator is adjusted 
to provide a 45 kHz output. The clock output consists of a sequence of 
consecutive binary one's and zero's to the PN sequence generator. The one 
state is +5. 0 Vdc and the zero state is ground level. 

PN Sequence Generator . The sequence generator is an ML, 20-stage 
shift-register with feedback from taps 17 and 20. 

Delay Circuitry . There is a delay circuitry (channel) which consists 
of a NAND gate matrix for coarse delay selection and clock vernier for fine 
delay control. 

Spectrum Shaping Filters . The shaping filters are matched. Butter- 
worth active lowpass filters. Filter insertion loss is zero dB, with zero 
dB ripple in the response band. 

Conditioning Amplifiers . Channels 1 and 2 conditioning amplifiers 
have gain and dc offset potentiometers (front panel) provided for control of 
the amplifier outputs. 

Difference Amplifier . The difference amplifier consists of three 
Fairchild 741 operational amplifiers. The amplifiers present a high imped- 
ance to the differential inputs, and has a 0. 01 ohm output impedance. The 
difference amplifier has a gain of 2; therefore, the output is 2 times the 
difference of the 2 inputs. 

Recorder Interface . Interface circuits are provided for the record- 
playback modes of operation. In the record mode, the outputs consist of a 
differentiated version of the clock and pseudorandom binary sequence. In 
playback mode, the interface accepts the differentiated signals from the 
recorder and decodes the clock and binary sequence signals. These signals 
must have an amplitude of at least 1.3 V at the PN test set input during 
recorder playback. 
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Figure 2. PN test set block diagram 
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Figure 3. 19 in. rack PN test set 


Electrical Specifications 


Inputs 


1. Analog response signal from test system. 

2. External clock-internal clock is 45 kHz. If other frequencies are 
desired, external clocks must be used. The unit has been tested with external 
clock frequencies up to 1 MHz. However, with the PN sequence generator used 
to ensure a Gaussian amplitude probability density distribution, the w, / co q 

should be 12 to 36 (allowable range) . 

Outputs 


1. Analog signal (filtered PN sequence) to system under test, vari- 
able 27 mV (p-p) to 5.5 V (p-p) with adjustable dc offset from -5.0 Vdc to 
+4.2 Vdc. The power density spectrum envelope of the analog signal is 
primarily that of the lowpass Butterworth shaping filter (the unit has select- 
able filter pairs with cutoff frequencies of 400 Hz, 500 Hz 1 kHz, and 2 kHz) . 
The amplitude probability density is near-Gaussian with an adjustable mean. 

2. Input/output average delay (back panel) vernier incremental 

commands to counter and front panel 21 position rotary switch 11. 5 ns to 

455 ps, readable to less than 0. 1 ps. 

3. PN Sequence — from ML PN sequence (20 states) when in record 
or playback modes (not normal test modes) . 

4. RMS error output signal — error output from difference amplifier 
provides interface with true rms (T-rms) voltmeter. 

Necessary Auxilliary Equipment . T-rms voltmeter, elapsed time 
counter (for vernier delay) . ” ' 

Power . 115 Vac, 60 Hz, single phase (less than 3 amps) powers 
internal 6 Vdc (1.5 A) and ±15 Vdc (0.2 A) power supplies. 
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IV. DATA FILTERS 


Active Data Filters 

The pseudonoise test set uses a matched active filter pair, one filter 
in the output signal channel and one in the reference channel. The filters 
are matched to minimize the true rms value of the difference in the response 
functions for a random sequence input. Five interchangeable plug-in filter 
sets were constructed. Four are Butterworth third-order lowpass filters 
with 3 dB cutoff frequencies at 400 Hz, 500 Hz, 1000 Hz, and 2000 Hz, and 
one is a first-order lowpass with 3 dB cutoff at 2000 Hz. The filter order and 
cutoff frequencies arc listed below: 

Filter Order (pole) Cutoff Frequency (Hz) 

3 400 

3 500 

3 1000 

3 2000 

1 2000 

General Design Theory . The filters were synthesized using the active 
network design procedure given in Computer-Aided Filter Design Manual, 
section 6.0 [ 1] . The computer programs provided in the manual were used 
to design the filters. The filters consist of active second-order (quadratic) 
lowpass networks. Odd order filters are realized by connecting a first-order 
lowpass network to the output of the last quadratic network. The quadratic 
network uses a voltage controlled source type of amplifier having a gain close 
to unity. In the PN test set, an operational amplifier at unity gain was used. 
The controlled source type of active network requires the minimum number of 
circuit components. 

The general voltage transfer function of a quadratic lowpass filter is, 



( 7 ) 



A quadratic (2 pole) lowpass network and a 1 pole lowpass network are 
connected in cascade to obtain the third-order filter configuration shown in 
Figure 4. This configuration was used to realize all the third order filters. 
The variable K is the gain of the amplifier, R 0 is the output resistance of 
the amplifier, and R is the load on the filter. 

Li 

The normalized transfer function of the quadratic section is 


Eg (s) _ Gs 2 + Ds + K 
^ As 2 + Bs + 1 


where 

B = RjCj (1 - K) + C 2 (R t + Rg) + R 0 C! , 

G = Ro 1^2 Cl c 2 , 

D = Rq Ci • 


(9) 

( 10 ) 

( 11 ) 

( 12 ) 


Here, R t includes the source resistance to the network and C 2 
includes the input capacitance of the amplifier. 

c, 



Figure 4. Third-order active lowpass filter. 
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The normalized transfer function of the first-order section is 


E 3 (s) = E 
E->(s) II s + 1 


( 13 ) 


where 


E 



R 3 + r l 


H 


^ 3 


(14) 


(15) 


For the relatively low cutoff frequencies involved and for K » 1 and 
R 0 « 0, the network transfer function of equation (8) reduces to. 


Ms) _ 1 

Ej(s) R t It, C,C 2 s 2 + C 2 (R t + Rg) s + 1 


(16) 


The synthesis equations for the networks described by transfer functions 
(16) and ( 13) sealed in frequency by cu are 
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where c o^ equals the 3 dB cutoff frequency in radians per second, R t R2 C t C 2 

and H are set equal to unity, and B = C 2 (Rj + R2) • The usual synthesis 
procedure for the quadratic network is to, first, select a value for C 2 and, 
then, choose C t such that Cj 2 4C 2 /B 2 . Then R t and Rg are calculated. 
For the matched filters, the more accurate synthesis equations of the computer 
program were used where K= 0.99999 and Ro=0.01 ohms. A thorough 
development of the exact design equations used is given in a thesis [2] which 
was the source for much of the active filter material contained in Computer- 
Aided Filter Design Manual [l] . 

The normalized pole positions for a third-order Butterworth filter 
are s= 0.5000 ± j0.8600 and s= -1.000. These poles give a lowpass 
transfer function of 


F(s) 


1 . 1 

s 2 + s + l s+1 


( 21 ) 


To frequency-scale the function to have its 3 dB cutoff at co^ , the 

complex operator s is divided by co^ . The transfer function of equation (21) 

scaled to have a 3 dB cutoff at co is, therefore, 

o 




o 


S 2 + CO S + CO 2 

o o 


CO 

o 


S + CO 

o 


( 22 ) 


Reference to equations ( 13) and ( 16) show that coefficients il and B 

are equal to unity and that the 3 dB cutoff frequency co is equal to the 

o 

natural frequency co^ of equation (7) . 

Matching and Tuning Procedure . An advantage of the active filter 
design procedure contained in Reference 1 is that tuning can be eliminated or 
at least held to a minimum. However, because of the extremely close 
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tolerances required to match the filters to less than 0. 1 percent rms error, 
a rigorous component matching and tuning procedure was necessary. The 
following procedure is used. 

i 

Two sets of approximately matched components are obtained for the 
values calculated by the computer program. Using one set of components, 
a filter is breadboarded and slightly tuned, if necessary, to have the specified 
cutoff frequency. The components are then removed from the breadboard 
and individually placed on a General Radio Type 1608-A Impedance Bridge. 

For each component, the bridge meter is nulled. Then, without disturbing 
the null setting controls, the component is removed and its approximately 
matched counterpart is placed on the bridge. Without adjusting the bridge 
controls, the component is then built up or adjusted until the meter reaches 
the same null. The same procedure is applied to all the components. 

A pair of closely matched operational amplifiers is selected from a 
large group of amplifiers. The amplifier parameters measured for matching 
purposes are the input capacitance and the output rise time. The filter pair 
is then constructed using the closely matched component pairs. At this point, 
the rms error for the third-order, 200-Hz filter was about 0.4 percent. This 
error is probably due to slight differences in the wiring layout, the amplifier 
parameters, and the resistive and capacitive components. 

To further reduce the rms error requires very fine tuning of the con- 
structed filter pair. Faeli component of the filter pair is individually tuned 
while observing the; effect on the rms error. First, each resistor is individ- 
ually paralleled with a 20 to 50 resistor and the effect noted. If the error 
is reduced, a parallel resistor is permanently added. Then each capacitor 
is parallel with a 1 to 5 pF capacitor and the effect is noted. This procedure 
is repeated several times until minimum rms error is achieved, usually below 
0.2 percent. At first the rms error drops only slightly for each component 
adjustment but begins to drop more rapidly as a more exact filter match is 
approached. Finally, each resistor is individually paralleled with a capacitor 
and each capacitor with a resistor, and the effect on the error is noted. This 
further tuning reduced the rms error of the third-order, 200-Hz filter from 
0. 1 percent to 0. 05 percent. 

An advantage of using active filters should be that the second- and first- 
order sections can be matched separately and then connected. However, this 
is not the case because of the extremely high degree of matching required. 
There is still a small amount of interaction between the section which must be 
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included during the tuning process. When the sections were matched separately 
and then connected, the combined sections were no better matched than before 
tuning. 

Design Results . The first matched filter set built was a third-order 
Butterworth with 3 dB cutoff at 2000 Ilz. The schematic for this filter pair 
is shown in Figure 5. The responses of the filter pair are matched to 0. 05 
percent true rms error. The total set error with the filter plugged in was 
0. 15 percent. The rest of the prototype test set introduces about 0. 1 percent 
rms error mainly due to radiation of digital switching spikes into the data 
and differences in output rise times of digital circuits in each channel. This 
error could be reduced by better separation and shielding of the analog output 
from the high speed digital circuitry. 


10 093pF 



Figure 5. Third-order, 2000-Hz matched filter pair. 
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Four other filter pairs were eonstructcd having the same configuration 
as shown in Figure 5. Of course, the one pole filter consisted of only R 3 
and C 3 . For these four filters, the fine tuning was accomplished by placing 
trimmer capacitors across all capacitors of one filter pair. By using this 
procedure, all capacitive tuning could be performed on only one filter as the 
capacitance could be easily increased or decreased. It was found that it was 
more difficult to reduce the rms error using the trimmers than with paralleling 
discrete components. 

The component values used for the other four matched filters along 
with their final rms errors are given below. 

For 400 Ilz cutoff, third-order filter: 


Il t - 78.596 KS2 
It, - 21.552 Ktt 
R :! = 27.491 IU2 
Filter rms error = 
Total test set error 


Ci = 25 525 pF 
C 2 - 5967.8 pF 
C 3 = 15 448 pF 
0. 12 percent 

= 0.228 percent (with filter plugged in) 


For 500 Hz cutoff, third -order filter: 


Rj - 67.115 K12 
1^2 = 17.067 KS2 
R 3 .= 20.9 KS2 
Filter rms error = 
Total test set error 


C t = 25 594 pF 
C 2 - 5776 pF 
C 3 = 15 880 pF 
0. 05 percent 
= 0.15 percent 


For 1 KHz cutoff, third-order filter: 


Ri .= 62.577 KS2 Cj = I 1 957 pF 
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[*2 If5. 95') M2 2001 pK 

Rg -- 20. 889 M2 C :t 7 997 pE 


Total rms error ~ 0.09 percent 
Total test set error ~ 0. 207 percent 


For 2000 Hz, first-order filter: 


R 3 = 20.054 Kfl 
Filter rms error = 
Total test set error 


C 3 = 3787.2 pF 
0. 08 percent 
= 0. 10 percent 


Sensitivity Analysis . Because of the difficult and time-consuming 

tuning required for construction of the four matched filter pairs, a sensitivity 

analysis was performed on the active quadratic networks. The damping ratio 

(O and the natural frequency (co ) of the network were taken as being repre- 

n 

sentative parameters of the shape of the response curve. Therefore, the 
sensitivity of t, and co were determined for a percentage change in each 

of the network components. The results indicate that the sensitivity of t, 
to changes in Rj or Rg is minimum when Rj = Rg, which occurs when 
Cg /C i = t, 2 . The analysis was performed for the simplified case where 
K = 1 and R 0 = 0. For designs where R 0 + 0, Ri must be somewhat larger 
than Rg. 

Equating the transfer function of the general quadratic function, 

equation (7) , with the quadratic network function, equation (16), where 

coefficient B ~ 2L, gives £, and c o as 

n 


£ 



• 1*1 + 

2 \l R, Rg 


CO 

n 


CO, 


^ Rj Rg C, Cg 


(23) 


(24) 
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where co equals lowpass 2 dB cutoff frequency. Normally R t • R 2 • C ( • C 2 
is made equal to unity and, therefore, = u^. 

The equations derived for the change in 4 and u)^ for percentage 

changes in components Rj, R 2 , Cj, C 2 and ratio C 2 /Ci are given below. 

The ratio of C 2 /Cj determines the separation in the values of Rj and R^ 
When Cj/C^ 4 2 , R t - IV 


34 


3R 


L 

Ri 


I &2 Hi ~ ^2 
Cl 4 Rj Ra 


a4 




f±L J .^2 % 

Cl 4 n/ R| Ra 


a?. 


ac 

Ci 


£ 

2 


84 

aCj 

Co 


£ 

2 


34 

9z 


4 

2 


8u) 

n 



OJ 

n 


3co 

n 

8R2 

*2 


0) 

n 

IT 


(25) 

(2(5) 

(27) 

(28) 
(29) 

( 20 ) 

(21) 
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9ca ca 



Ci 


The previous equations show that the change in ca for a percentage 
change in R t or Itg is minimum when Rj equals R^. The change in c a ■ 
for a percentage change in either Cj or C 2 is independent of the component 
values. However, the change in t, for an absolute change in Cj or C 2 
is inversely proportional to the value of Cj or C 2 . For example, 

_9£ 

9C 2 ' 2C 2 


Hence, C 2 should be made as large as possible for ease in making a small 
change in £,. 

The change in ca for a percentage change in any one of the components 

is independent of any of the components' values. However, the change in ca^ 

for a change in value of any one component is inversely proportional to that 
one component. For example. 


9ca 

n 

9R7 


ca 

n 


2R t 
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Since the filler is designed by first selecting C 2 and then selecting 
the ratio CjC^ & l/£ 2 , it appears that the best design for easiest tuning 
would be to select C 2 as large as possible. Then choose C t equal to 
C 2 A 2 which makes Rj equal 1^. For the four filter pairs constructed, 
the value of Rj was much larger than R^ The tuning might have been 
somewhat easier if R t and R^ were approximately equal. 


Passive Data Filters 

The application of passive data filters for the pseudonoise test set has 
been inhibited because of the difficulty in tuning these type filters. Unlike 
the active configuration, each filter pole position in a passive filter depends 
upon the interrelationship of all the network components. This complicates 
the tuning procedure because the relationship between critical filter network 
characteristics and the network components can be quite involved. For the 
purpose of designing passive filters for pseudonoise system testing, an 
improved filter tuning procedure has been developed. With this procedure, 
effort in tuning a particular filter element is in proportion to that element’s 
effect upon the critical network parameter. The application of passive data 
filters for PN testing might be desirable if the tuning procedure developed for 
the passive design results in a very accurately tuned filter. 

The tuning procedure developed uses the rms value of the difference 
in the response functions of a filter and its theoretical perfect reference as 
the critical parameter. This rms difference is calculated as a function of 
the variation of each filter element from its synthesized value. Assuming 
that the rms difference function is accepted as the critical parameter in 
filter design, then those elements whose variation most affect this parameter 
should receive proportionally more effort to assure their accuracy. 

Mathematical Model . Figure 6 is a block diagram of the basic model. 
The uniform spectral generator drives each filter with a waveform that has a 
constant power spectral density. Each filter weights and delays the spectral 
components, and the normalized rms value of the difference waveform is 
the critical parameter. 
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Figure (i. Basie model block diagram. 


The normalized mean—squarc error can be expressed as 


CO 

C 2 

J {Fi (cu) I COS [cut - 0} {«)] - I F 2 (cu) | cos [cut - 0 2 (cu) ] } dcu 


CO 

°° J [ F t (cu) I 2 cos 2 [cut + 0 t (cu)] dcu 

— CO 


(.‘55) 


In expression (35), the over-bar indicates time average, and other terms are 
defined as follows: 

Fj (cu) transfer function of the first (reference) filter. 

F 2 (cu) transfer function of the second (modified) filter. 


Define: 


Fl ^ ' D t (cu) Re [Dj (cu)] + jlm [D t (cu) ] I F 1 ( w ) I 6 * ^ ^ ( 36 ) 
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F, (co) 


D, 


(co) Rc f I) 2 (co)1 + jlm fD 2 (co) 1 


= I F, (co) I G 


jOo (co) 


(37) 


where 


Oj (co) = - tan 


0 2 (co) = - tan 


| t p i (u) | \ 

\ Re [D, (co)] j 

j Im [D, (co)] \ 

( Re [D 2 (co) ] j 


(38) 


(39) 


and D, ( co) and I) 2 (co) are complex polynomials in co . Expanding the 
expression taking the time average yields 


i'<- 


(J 


(T 


2 I - orj 


J I Fj (co) I 2 dto - 2 | I Ej (co) I I F 2 (co) I cos f0 1 (co) - 0 2 (co) I dco 


■l / I I 1 ', (co) I 2 dco 



I I F t (co) I 2 dco 


(4 0) 


or 


2 ( 00 ® 

“ \ / t F t (co) | 2 dco - 2 J I Fj (co) | |F 2 (co)| cos ]0j (co)] cos T0 2 (co)] dco 

V ( lc0 

CO 

-2 f |F](co)| i F 2 (co)I sinlO^co)] sin[0 2 (co)] dco 


-i f | F 2 (co) | 2 dco 



f I F { ( co ) ! 2 dco 


(41) 
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Substituting into equation (4l), the relations 


COS [0J 

(co)] 

- Re [ D, (co) 1 | 

I' l (co) 1 

(42) 

sin [0 l 

(co)] 

= - Im [D, {co) ] 

1 F ! (co) I 

(43) 

cos [0 2 

(co)] 

= Rc [D 2 (co)] | 

F 2 (co ) | 

(44) 

sin f 0 2 

(co)] 

= Im 1 D 2 (co) 1 | 

F 2 (co) | 

(45) 


yield 


= < J | F^co) | 2 do) - 2 j | F t ( co ) | 2 | F 2 (co) I 2 Re[Dj(co)] Re[D 2 (co) ] dco 


- 2 f | F 4 ( co ) | 2 | F 2 ( co )| Irnf D^u))] • Im[D 2 (co) ] dco 


-i f | F 2 (co) | 2 dco 


/ | Fj(co) | 2 dco 


= (1 2 - 2I 3 - 21, 4 l 2 )/l t . 


(46) 


Equation (46) may be used to perform an element sensitivity analysis 
for networks whose transfer functions can be expressed as 


F(w) = 


1 

Ij(w) 


(47) 


where w is some function of co . The use of w allows the analysis of 
bandpass, lowpass, and band-stop filter configurations by transforming to a 
lowpass equivalent form, equation (47). 



Fil ter Specialized to Cauer Configurations . Figure 7 contains block 
diagrams of the Cauer singly terminated networks for odd and even orders 
respectively. Figure 8 shows a section of the ladder network where the 
elements have been specialized to a lowpass configuration. 


The dissipative and reactive components in each element may be lumped 
together by defining 


jeo* = jeo + d 


(48) 


where 


d = l/Q , (40) 

and d is called the element dissipation factor, Q is the element quality 
factor, and ca* is called the dissipative frequency term. 


F(m) F(m-2I F(m-4) 



F(m-1) 


a. Odd order. 


F(m) F(m-2) F(m-4) 



F(m-1) 


b. Fven order. 

Figure 7. Cauer configuration. 


c- 
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L j+2 R j+2 



Figure 8. Section of a ladder network, specialized 
to a lowpass configuration. 


Realization of Chebyshev and Butte rworth Filters . The form of the 
Butterworth transfer function is 


F(w) 

\/ 1 -I' LO 


( 50 ) 


where N is the filter order. The form of the Chebyshev transfer function is 


F(a>) - ■ — 

n/ 1 + e 2 r I’ N 2 (w) 


( 51 ) 


where e is a real constant, e < 1, and T n (ca) is the nth order Chebyshev 
polynomial defined as 


t n M ' 


cos 


(n 


cos 


- 1 


Z) 


( 52 ) 


where Z is a real variable, Z 1 . 


Reference 1 lists a group of programs that synthesize filter elements 
for Chebyshev and Butterworth filters. These programs perform the following- 
functions: 

1. Using an expansion of the Chebyshev polynomial calculate the poles 
of F(a>) . 

2. Predistort the pole plot to accommodate dissipative components 
in the filter elements. 

3. Perform a synthetic, division to calculate ladder (Cauer configura- 
tion) element values. 

The actual synthesis procedure involves the following steps: 

The transfer omittance of a terminated ladder network is 


h. - _ 

Vi 


Yl2 


.. Y 12. G 2... 

G 2 + Y 22 


(53) 


where G 2 is the admittance of the termination. If G 2 is normalized to 1, 
this equation takes the form 


- Y 


- Y 


12 


12 . 


• + y 22 


(34) 


The transfer function of the network is 


Gl2 


— Y tt 
G 2 + Y 22 


(55) 


and in the case of normalized load impedance, it is seen that 


G 


12 


- Y r , 

1 + Y 22 


( 5(5 ) 
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The transfer function can be written in the form 



P (s) 

qi (s) + q 2 (s) 


(57) 


where (s) is a polynomial composed of the even powers of the polynomial 
in the denominator of the transfer function, and q 2 (s) is the polynomial 
composed of the odd powers. Also, equation (57) can be rewritten as 


P (s) 

p ( h) = q (s) 
qi (s) + q 2 (s) 1 , q? (s) 

qj (s) 


(58) 


and 


Y 22 


tfe («) 


(59) 


If the choice of (s) and q 2 (s) was not made in the manner given, 
the network would not be realizable as a grounded structure. 

Figure 9 is a block diagram of the transfer function calculation and 
filter element modification portion of the program. The first block represents 
the transfer function calculation from the Chebyshev polynomial expansion. 

The second block includes the isolation of the output admittance Y 22 by 
assuming a normalized resistive termination. 

The parallel paths include the calculation of the transfer function for 
both the reference filter and the modified filter. This procedure may appear 
redundant because the transfer function for the reference filter is available 
in block number 1 of Figure 9. The procedure given was selected to minimize 
the effect of round-off error in the numerical calculations. 

The transfer functions are calculated from the modified and unmodified 
network elements by solving a matrix equation of the form 


KT - K 


(60) 
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Figure 9. Transfer function calculation and filter element 
modification block diagram. 











where 


I = 


( — + 1 1 order vector, 

A 2 / 


m even 


order vector, m odd 




i 1 1 order vector, 

V 2 / 

N- 1 ) 


m even 


order vector, m odd 


where E< =£ 0, E. - 0, \ + \ , 
i 


+ 1 ) order matrix, m even 


K = 


(?•■) 

N-0 


order matrix, m odd 


T 

The K matrix is symmetrical about the diagonal, K = K, and 

k - 0, i - 1 > i > j + 1, where k is the iith element of the K matrix, 
ij i.i 

Example for a First-Order Filter . A first-order filter in the config- 
uration synthesized by the method of Reference 2 is shown in Figure 10. 

The dissipative frequency term 


co* = co - jd 
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Figure 10. Cauer configuration for a first-order filter. 

was defined in equation (48) . As an example of the type of calculations per- 
formed by the sensitivity programs, a first-order example is given here. 
The transfer function of the filter is 


G i2 




(j£*>) 


L 



( 01 ) 


where is the dissipative component of the network element, and Fj is 
the reactive component. 

To simplify the calculations for illustrative purposes, assume for the 

reference filter in Figure 10 Fj = 1, Z = 1, and dj = 0 . 

JL 

From Reference 8, 


a 


1,1 



( 02 ) 


a l, 2 


F 1 (I 1 + 




1 


9 


(o:i) 
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(<M) 


a 2, 1 



a 


2,2 


F t d, + Z, 


(65) 


The terms with bars indicate parameters from the F 2 filter in Figure 10. 
Referring to equation (46) , the integrals can be evaluated as follows. 

Evaluation of L : 


D l ” a l,2 “ 1 
D 2 = a lj2 = 1 

Ki = 1 

h = (gi/2l>i D 2 ) ; 1/2 


Evaluation of I 2 : 


Fid, + Z 1 


[)h = a. 


2,2 


Dn = a. 


2,1 Z 


gl = 1 


I 2 = 


Jl 


( 5 l) ! 


2Dl ° 2 2F,(F 1 d, + Z L ) 


( 66 ) 


(67) 
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Evaluation of I 3 : 


F 1 d t + Z. 


D l ~ a l , 2 a 2 , 2 


D 2 - 3^2 a 2, 1 + a 1 , 1 a 2,; 


- F< d« + Z 

V t 1 1 L 


Z 


D 3 = 


a t 1 a 2 ! - 

9 9 



Z 


L 


gl = 0 


§2 - a l,2 a 2,2 _ 


Ftdj + Z, 


L 


I = z-l Llh V ikJh 

J 2D t D 2 I) :t 


2 



Z 


L 


1 F i d t + Z J 


( 68 ) 


Evaluation of I 4 : 


F 1 dj + Z 


Dj - a^ 2 a 2j2 


F F i + Z L 

D 2 ~ a l s 2 a 2 1 + a l, 1 a 2 2 ~ Zr" + Z 

Z L Z L 


d 3 


‘1,1 a 2,l 



z 


L 
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gi - - a ( j aj ) 


Ij 


g 2 = 0 


I = - Si D < + fe Da 
4 2D, Dj D, 


2(F 


j + F j dj + Z . 


) 


(GO) 


From equation (GO) , the normalized mean square error is 


4Z. 


t = 1 - 


F j + F -j cl j + Z . 


+ 


Z L 


(Ft) (F, dj + zj 


(70) 


where t is equal to a 2 /a 0 2 . Asa chock of equation (70) note that if 

Z = 1, F t = 1, and dj = 0, then c - 0 . This is a perfect match with the 
-L 

reference filter. 

The three parameters that can vary in the first-order filter are the 

ig: 

1. dj, the dissipative component. 


following: 


1. 

d t , 

2. 

Z L’ 

3. 

Fi, 

Case 1: 

e = 

1 - 


4Z . 


Ij 


1 + Z 


+ z 


Ij 


(71) 


Ij 


The dependance of e upon Z for Case 1 is illustrated in Figure 11. 

J-/ 
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ERROR 



Figure 11. t us a function of Z . 

L 


Case 2: F, / f, 7, x - 1, dj - 0 


F, + 1 

The dependancc of t upon 
Case 3: dj ± 0, Fj = 


2 + dj 

The dependance of e upon 


Fj for Case 2 is illustrated in Figure 12. 

1, Z L .l 

1 

1 + dj 

d t for Case 3 is illustrated in Figure 13. 


(72) 


(73) 
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ERROR (e) ERROR 



NETWORK COMPONENT (F^) 

Figure 12. e as a function of Fj. 



The variations of e for the three cases given include a wide range of 
variations of the network elements. Normally, relatively close adjustment 
of each parameter can be achieved with a reasonable tuning procedure. More 
useful information in this case would be a parameter reflecting the sensitivity 
of e to each parameter near its specified value. A useful parameter for 
this purpose is the partial derivative of e with respect to the network 
element of interest. For the first-order filter, these include: 


9e 


8Z . 


Fi 



0( dj - 0 

8F t Z - 1 


(■•SFl 1 I) (Fj - 1) 
Fi 2 (F t + l ) 2 


Be 

Bdf 


F t -• 1 di (3d t + 4) 

Z = 1 (d t + 2) 2 (dj + l) 2 


(74) 


(75) 


(76) 


If the values of the network parameters are near their design value, 

Fj = l + A F [ , d t - Ad 1} and Z 1 -i Z for the example given, then 

L 


at 


az. 


L 


AZ 


L 


Fi 1 


1 i AZ 


I, 


(77) 


8e 

9Fj 


dj - o 

Z - 1 

L 


AF, 


I i 3AF 


(78) 


9e_ 

8d 


Fi = 1 
'L 


Ad, 


Z_ = 1 1 + 3Ad t 


( 79 ) 


where AZ, AF l9 and Ad t are small. 
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The tuning procedure developed should improve the performance of 
passive filters for applications where accurately matched transfer functions 
are required. The transfer function match requirements upon the data filters 
in pseudonoise system testing should be within the capabilities of passive 
filters if this tuning procedure is applied to the filter manufacture. 

V. PREPARATION AND OPERATION OF PN TEST SET 

Preparation 

This section discusses the PN test set installation, controls, and 
indicators. 

Installation . The PN test set chassis was designed for mounting in a 
standard equipment rack. A power cable is provided for connection to a 
115/125 Vac power source. All system interconnects are on the front panel. 
Two BNC connectors are mounted on the back panel to facilitate use of an 
electronic counter to measure exact vernier delay. 

Controls and Indicators . Controls and indicators on the PN test set 
front panel (Fig. 14) are shown in Table 1. 


Operation 

This section describes normal use of the PN test set to evaluate 
characteristics of a telemetry-communication system in both the real-time 
and record-playback modes of operation. 

Real-Time rms Error Measurement . 

Measurement. Measurement of real-time rms error due to waveform 
distortion in a communication system can be accomplished by the following- 
procedure. 

1. Connect the test setup as shown in Figure 15. 

2. Place the test set mode switch in the real time position. 


:!7 
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Figure 14. PN test set controls. 



TABLE 1. PN TEST SET CONTROLS 


Key 

Control 

Type 

Purpose 

1 

POWER 

Pushbutton switch 

Applies 115/125 Vac power 
to the test set 

2 

TEST MODE 

3 -position rotary 
switch 

Selects mode of operation 

3 

CLOCK 

Two-position 
tousle switch 

Selects internal or external 
clock 

4 

GAIN CH 1 

10-turn potenti- 
ometer 

Controls gain of the refer- 
ence signal during REAL 
TIME and PLAYBACK 
operation; Controls gain of 
data signal during RECORD 
operation 

5 

GAIN CH 2 

10-turn potenti- 
ometer 

Controls gain of the data 
signal during REAL TIME 
operation 

6 

DC OFFSET 

1-turn potentiom- 
eter 

Not used, control is now on 
card 6 

\ 

7 

DELAY INCRE- 
MENTAL 

21-position rotary 
switch 

Selects CH 1 delay in 22.2 
fis steps 

8 

DELAY VERNIER 
IN -OUT 

2 -posit ion toggle 

Switches VERNIER DE LA Y 
circuit into total delay 
circuit 

9 

DELAY CONTROL 

1-turn potentiom- 
eter 

Controls CH 1 vernier 
delay from -11. 1 ns to 
+11. 1 gs 

10 

RMS VOLTMETER 
ERROR X2 

BNC Connector 

Provides connection to an 
rms voltmeter for meas- 
urement of X2 error 












TABLE 1. (Concluded) 


Key 

— 

Control 

Type 

Purpose 

11 

FROM SYSTEM 

BNC Connector 

Input from system under 
test during REAL TIME 
and PLAYBACK operation 

12 

TO SYSTEM 

BNC Connector 

Data output to system 
under test during REAL 
TIME and RECORD oper- 
ation 

1 2 

PN SEQUENCE 

BNC Connector 

Output/input for PN 
recording track during 
RECORD and PLAYBACK 
operation 

14 

EXT CLOCK 

BNC Connector 

Output jack for clock signal 
recording during RECORD 
and clock input during 
PLAYBACK 

Back Panel Controls 


EXT CLOCK 

BNC 

Input jack for external 
clock during REAL time 
and RECORD 


START DELAY 

BNC 

Output jack for reference 
clock 


STOP DELAY 

BNC 

Output jack for delayed 
clock 
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Figure 15. Test setup for real-time rms error measurement. 

3. Adjust channel 2 gain control to the correct system dynamic 

range. 

4. Adjust the front panel channel 1 gain control, incremental delay 
and vernier control to obtain the minimum reading on the true-reading rms 
voltmeter. 


NOTE 


This cancels amplification and delay differences, 
and the rms voltmeter indication is 2 times the result- 
ing distortion. 

Theory. Theory of this measurement concept is based on the definition 
of rms error: 


e 

rms 


[e 2 (t)] 1 / 2 


(HO) 


True-readings rms voltage measuring devices give an indication proportional 
to an average power absorption over a period T., the instrument time con- 
stant. Specifically, 1 


e (t) = 
rms 



f e 1 (t) t“ k(t - T) dr 

0 


( 81 ) 
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If the approximation is made that 


- k(t - T) = r, t 


T < T. 


e k (t ‘ t) - 0 , t - r > T. , 


then 


e 

rms 



[g2 (t)1 T., 


)l/2 


(82) 


Figure 16 is a block diagram of the system internal to the PN test set. 
For the case shown in the figure, the true-reading rms volmeter indicates 


e rmg (*) = 1 1 e l (t) - <=2 (t - t ) ] x 2 | / ( 83 ) 



Figure 16. Error detection section of the PN test set. 

I 
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or 


e ™J( fc ) + t « 2 2 ( 1 - T ) 1 T -2[e i (t)e 2 (t -t)}^ . (84) 


rms 


The PN test set output is a twofold amplification of the rms error; therefore, 

the rms meter reading equals 2e 2 . 

rms 

Record-Playback rms Error Measurement . Measurement of rms 
error due to distortion in a record-playback situation can be accomplished 
with the PN test set. Data from the system under test and the pseudonoise 
sequence is recorded and then played back to obtain the original data and 
measure the rms error. 

Record Mode . The procedure for operating the PN test set in the 
record mode consists of the following steps: 

1. Connect the test setup shown in Figure 17. 

2. Place the PN test set controls in the following positions: 

a. Mode switch to record. 

b. Vernier delay switch to out. 

c. Incremental delay switch to position 1. 

3. Adjust the Ch 1 gain potentiometer to match the recorder interface 
requirements. 

4. Record data as required. 



J, TO SYSTEM CHI 

1 

PN 


TEST 

J 4 PN SEQUENCE 

RECORDER 



Z (TRACKS) 

SET 

Jg CLOCK 

3 


Figure 17. Test setup for recording. 
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Playback Mode . The procedure for operating the PN test set in the 
playback mode consists of the following steps: 

1. Connect the test setup shown in Figure 18. 

2. Place the mode switch in the playback position. 

3. Activate the playback on the recorder. 

4. Adjust the following PN test set controls to obtain the minimum 
rms error indication: 

a. Ch 1 gain. 

b. Incremental delay. 

c. Vernier delay. 


NOTE 


The rms voltmeter will read 2 times the rms error 
due to distortion in the record-playback process. 

Time Delay Measurement . Average time delay in a telemetry- 
communication system may be computed by the following procedure. 



Figure 18. Test setup for playback. 
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1. Connect the test setup shown in Figure 19. 

2. Using the procedure described in the section on real-time rms 
error measurement, measure the rms error. 

3. Using a counter or an oscilloscope, measure the vernier time 
delay (in Ch l) required to give this minimum rms error indication on the 
rms voltmeter. 

The distance between negative-going edges of the two pulse-trains 
(start output and stop output) is the vernier time delay (d) . 

4. Compute the total time delay as 

D t = (k - 1/2) T + d 

where k is the setting on the incremental delay switch and T is the clock 
period (22.2 /js for internal clock) . 

The setting of the incremental delay switch gives the additional delay 
in Ch 1, in steps of T (T equals 22.2 n s for internal clock operation) . 

When the vernier delay switch is the out position, is equal to (k)T 
seconds. 



Figure 19. Test setup for time delay measurement. 
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The measurements in the three previous subsections, the Ch 1 and 
Ch 2 dc offset pbtentiometers, mounted on card 6, need to be adjusted to 
obtain minimum rms error only when a dc reading true rms voltmeter is 
used. Most true rms voltmeters are ac reading only and do not include dc. 
The dc offset potentiometers are used to set the dc level of the data and 
reference signals. 

Correlation Function Determination . Determination of the correlation 
function of a communication system may be accomplished by the following 
computation. 

Since the correlation function for periodic waveforms is, by definition 


T 0 /2 

1 c 

Re 2 e t ( t ) =7 J e i M e 2 ^ “ r ) dT 

0 -T 0 /2 


(85) 


where T fl is the basic period of the waveforms. Assuming that 


T. = 

l 



then 


e „ 2 = t e i 2 (t) ] + [e 2 2 (t - r)] - 2Re 2 e t (r) , (86) 

rms i Q x 0 


or 


Re 2 e t (t) = 1/2 (e 2 + e 2 - e 2 \ . (87) 

y lrms 2rms rms J 

From the preceding equation, the correlation maximum occurs at the minimum 

rms error and by measuring system input, system output, and minimum rms 

error, the correlation of the system input and output is obtained. Remember 

that the PN test set measures 2e , and the value to be substituted into 

equation (87) is e 2 . rms 

rms 
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Impulse Response Determination . Determination of a telemetry 
communication system impulse response may be accomplished by the follow- 
ing procedure. 

1. Determine the system correlation function in accordance with 
the instructions in this section. 

2. Measure the system time delay (t . ) in accordance with the 

min 

instructions in this section. 

3. For values of delay different from , measure Re^ e t (t) 
and obtain impulse response from equation (91). 

Theory. Determine impulse response for the system with input, 
e (t) , and output e 0 (t) . 

1. The system correlation function is < 


CO 

Re. e 0 (t) = f e. (t) e 0 (t - t) dt (88) 

- CO 


also 


CO 

Re. e 0 (t) = J Re. e. (t) h (t + r) dt (89) 

-CO 


Where h (t + t) is Green's function. 

2. If the input is white noise, or looks like white noise to the system 
under test, then 


Re. e. (t) = kS (0) , (90) 
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and 


Re. e 0 (r) = kh (r) 


(91) 


Therefore, assuming the system under test is linear, the cross-correlation 
function is proportional to the system impulse response. 

VI. ACCURACY ASSURANCE OF THE PN TEST SET 
Interrelationship between Time Delay and Gain 

In an analysis of communication system performance [4] , the mean- 
square error between the input and output of a telemetry system has been 
expressed as 


a 2 (t, G) = lim 

T— *-co 


1 

2T 


JL 

/ [e. (t - t) - Ge 0 (t) - Gn 0 (t)] 2 dt (92) 
-T 1 


where e (t - r) is the input signal delayed in time, t is the time delay, 

G is the system gain, e 0 (t) is the output signal exclusive of noise, and 
n o (t) is the system noise. The waveform distortion is defined as the 
minimum of equation (92) which is obtained by proper selection of the vari- 
ables r and G. The values of r and G that result in minimum error arc 
defined as the average time delay, t 0 , and the average system gain, G 0 . An 
important property to be considered in the measurement of waveform distor- 
tion, both analytically and experimentally, is that the parameter t 0 can be 
determined independently of G () . This property is illustrated by replacing 
equation (92) by the equivalent frequency domain expression (ignoring the 
noise term) given by 


a 2 (t, G) '= f s(e>){l - 2GA(w) cos [0(cu) + cut] + G 2 A 2 (co)}da> , (93) 

0 
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where s(co) is the input signal power density spectrum, A(o>) is the signal 
amplitude response, and 0(a>) is the system phase response. The values of 
r and G that minimize equation (93) are determined from the equations 


9a 2 (t, G) _ 
9 r 


(94) 


9cr 2 (r , G) 
9 G 


(95) 


Performing the integrations of equation (93) gives 


CO 

0 = f o)s(cj) A(co) sin [ 0 (oj ) + car] do) , (96) 

0 


and 


G 


f s(co) A(oj) sin [0(w) + ut] do) 

_0 

00 

f s(a>) A 2 (w) dw 
0 


(97) 


Hence, as noted from equation (96), solutions for r may be obtained 
independently of G. If more than one solution is determined from equations 
(96) and (97) , the desired solution is the one which results in minimum 
error. 

The independence between r and G is illustrated in Figure 20 
where rms error is plotted as a function of time delay introduced into a band- 

limited random signal by an RC lowpass filter. The error js found to be a 

minimum at t 0 =. 0.046 msec independent of the selection of G. The curve 
G = 1. 124 represents the error for the average gain found in equation (97) , 
while the minimum of the curve defines both the average time delay and the 
waveform distortion. 
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Figure 20, RMS error as a function of time delay introduced 
into a band-limited random signal by an RC filter. 

The fact that the average time delay can be determined independently 
of the system gain suggests that, in a laboratory test unit for measuring 
waveform distortion, one should first seek to minimize the error by adjusting 
the delay and then determine the absolute error minimum through adjustment 
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of the gain. The PN test set, with indepenently designed gain and delay 
capability, readily lends itself to this procedure. It should be mentioned 
that, in using this equipment, some adjustment back and forth between gain 
and time delay may be needed to achieve the desired measuring sensitivity. 

This would be the case where the initial gain setting is far from the desired 
setting, causing the error versus delay curve to be almost flat. 

Experimental Accuracy Results 

The total measuring accuracy of the PN test set has been 
examined in two ways. First, the noise floor of the unit was determined 
by simply connecting a coaxial cable between the system input and output 
terminals. It was found that the unit could be nulled to such a point as to 
insure less than 0.2 percent rms error with data signal levels ranging from 
1.5 V p-p to 5 V p-p. Secondly, the distortion resulting from a more real- 
istic blackbox configuration was determined and the measured values were 
compared with calculated results. The following information describes this 
second method for checking the accuracy of the PN test set. 

The accuracy of the PN test set was checked against calculated data 
for several lowpass filters serving as blackboxes or test units. In general, 
the results substantiate the argument that the PN test set ( in its use of 
determination of end-to-end rms errors) is functioning properly. Although 
all of the experimental values taken did not match the calculated values exactly, 
the variations were small, and explained. 

There were several test filters (blackboxes with known characteristics) 
investigated with the PN test set. The basic setup used is the one shown in the 
block diagram in Figure 2. 

In all analytical work, the calculated data were shaped from pure white 
noise which had a uniform density spectrum and a Gaussian amplitude prob- 
ability density. All calculations were based on unmodified gains (i.e. , if the 
unit under test changed the amplitude overall, this was not adjusted for 
minimum error) . 


Two characteristics of the PN sequence-generated data (as opposed 
to filtered white noise) should be explained. 

First, consider Figure 21a. The calculated data were shaped from 

pure white noise and were assumed flat to co , whereas the measured data 

o 

(Fig. 21b) is 3 dB down at the One would expect less error in the 

measured than in calculated as lower amplitudes were in the w region. 
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WHITE NOISE SPECTRUM 






Second, the calculated data were shaped from pure white noise as in 
Figure 21a. Thus, the filtered data are assumed flat up to c and is then 

attenuated at 6 m dB/octave (m is the number of poles in the filter) above 
cutoff. The PN sequence-generated data, from which the empirical data 
results were taken, only approximate the properties of white noise. In con- 
trast to pure random noise, the PN sequence is actually periodic (refer to 
Figure 21b) and its spectrum nulls at integral multiples of co. , the clock 

frequency of the PN sequence generator. 

As illustrated in Figure 21, the difference in spectral characteristics 
can be significant for ratios of co /co = 20 for a shaping filter with low stop- 

K O 

band attenuation. For example, with co. /co = 20 and a filter attenuation rate 

k o 

of 6 dB/octave, the data in the vicinity of the first PN sequence null, co, , :is 

k ; 

only attenuated about 27 dB. Thus, for the calculated errors, data compo- 
nents exist at co although attenuated 27 dB, while for the measured errors 

K 

no data components exist at this frequency. Consequently, there is a differ- 
ence between calculated and measured data in the vicinity of the PN sequence 
nulls as well as the difference noted near the filter cutoff frequency due to 
the characteristics of the approximated spectrum. 

The way the data were taken was simply to vary co^ of the test unit 

filter and take the rms error readings with the test setup as is shown in 
block diagram in Figure 2. High order (sixth) filters were used as well as 
lower (first) order filters in testing. Figure 22 shows the very close corre- 
lation of calculated and measured results for all filters when a third-order 
Butterworth data shaping filter pair was used. As one would expect, better 
correlation can be obtained by using higher order data shaping filters and 
blackbox or test filters (for the reasons explained above) . 


VII. TECHNICAL PROBLEMS ASSOCIATED WITH THE 
REALIZATION OF THE PN TEST SET 

Two important technical problems surfaced in the design, construction, 
and evaluation of the PN test set. They included the extreme measurement 
accuracies required to tune the data filters to a mismatch of 0. 1 percent or 
less, and a problem with skewing of the amplitude distribution of the pseudo- 
analog noise signal for large ratios of sequence clock frequency to data filter 
cutoff frequency. 
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Figure 22. Plots for accuracy tests of PN test set. 


Filter Match Problem 

In the case of the filters, an initial passive design lowpass filter pair 
was matched to 0.4 percent rms error. An active design was developed and 
matched to better than 0. 1 percent. Section IV gives the details of the active 
data filter set design and analysis. 

For the purpose of developing an improved design procedure for 
matched data filters for pseudonoise testing, a filter component sensitivity 
analysis was performed. The analysis covered passive Butterworth and 
Chebyshev approximation type filters. The result was a mathematical pro- 
cedure for isolating the error sensitivity of a component on comparison basis 
and for standardizing the tuning procedure. Section IV gives the results of 
the analysis and an example of its application. 

Amplitude Distribution Skewing Problem 

For the particular maximal length digital sequence generator designed 
into the PN test set, the amplitude distribution skewing problem was observed 
to occur for ratios of clock to data filter cutoff of 50 and above. Ideally, the 
distribution should approximate a normal curve as illustrated in Figure 23. 
This was found to be the case for the lower ratios of clock to filter cutoff 
frequency; however, for ratios of 50 and above the distribution was skewed as 
shown in F igure 24 . 
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Figure 24. Skewed density function for large ratios of sequence 
clock to data cutoff frequency. 


An investigation of the causes of this skewing was undertaken in an 
effort to correct the problem. A literature search revealed that R. P. Gilson 
observed the phenomena in question and reported the results of an experimental 
program concerning filtered maximum length sequences [5]. Gilson’s report 
confirmed that the skewing problem was basic to filtered ML sequences, and 
not just a property of the particular ML generator in the PN test set. James 
H. Lindholm reported the results of an analysis of the statistics of M-tuples 
from ML sequences [6] . If the impulse response of the data filter is approx- 
imated as a weighted sum of the last M-bits from the sequence generator, 
Lindholm’ s results can be applied to the problem in question. Lindholm 
showed that the third moment of the distribution of weights of M-tuples from 
a maximum length sequence is related to the number of trinomials of order 
M-l or less that contain the sequence characteristic equation as a factor. 

In the case of the PN test set, the M-tuple length corresponds to the ratio 
of sequence clock to data filter cutoff frequency. The ML sequence generator 
in the PN test set is a 20-stage unit with feedback from the 20th and 17th 
stages. The sequence characteristic polynomial is 


X 20 + X 3 + 1 


(98) 



and this polynomial apparently is a factor of many trinomials of order 50 or 
less. An investigation was made to determine if there are any sequence 
generators with approximately 20 stages that would have a characteristic 
polynomial that is a factor of few trinomials of order 50 or less. It was 
determined that the 23 -stage sequence generator with feedback from stages 
23, 17, 11, and 5 has such a characteristic ploynomial. A computer simula- 
tion of the filtered sequence, with a filter impulse response period of 60 clock 
periods, was performed. Figure 25 is the result of this simulation. The 
solid curve is the approximated normal distribution and the • marks give the 
weights from the simulation. The approximation is satisfactory for the pur- 
poses of pseudonoise testing, and this sequence generator is planned for 
future models of the PN test set. 



Figure 25. Normal probability and observed weights. 
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VIM. CONCLUSIONS 


Although rigorous specification of the characteristics of filtered PN 
sequence data is apparently quite involved, the PN test set shows that a use- 
ful, pragmatic unit is available to provide rapid end-to-end rms error testing 
of communication systems. This unit uses filtered pseudorandom data and 
an analog error measuring approach to provide near-Gaussian, wideband data 
as a test signal. The test signal provided by the test set exercises the entire 
designed spectrum of the communication system under test. 

This paper explains the concepts of system testing of both wideband 
random and simple, deterministic data as well as the problems involved in 
utilizing these digital and analog measuring techniques. The theory of PN 
testing is outlined and the description of the PN test set effectively shows 
how the theory is implemented. Its development has led to improved methods 
for fabricating and accurately matching filter pairs. Analytical and empirical 
tests were conducted to show that this unit can accurately measure system 
performance characteristics. 

Two basic problems were noted during the project investigation: the 
need of a pragmatic method to overcome the extreme measurement accuracies 
needed to closely match the data shaping filter pairs, and the deviation of the 
data amplitude probability density distribution from the desired near-Gaussian 
distribution. Filter sensitivity analyses were used to overcome the filter 
fabrication problem and, while a specific PN sequence with a limited ratio of 
filter cutoff to clock frequency solved the immediate problem of a non- 
Gaussian distribution for the PN test set, additional investigations are being 
performed for more advanced data simulations for communication system 
testing. 


58 



APPENDIX A 


EXPERIMENTAL METHODS (ANALOG AND DIGITAL) FOR 
MEASURING WAVEFORM DISTORTION 


The basic intent of this appendix is to report two ways to perform 
end-to-end rms error testing on a system. Both schemes employ the idea 
of feeding a known signal into a system and comparing the system output with 
the known input to determine the degree of distortion. The question is how to 
data -load a system to get the most accurate test. The answer to this question 
must remain subjective until the specific system and its intended input data 
characteristics are known. However, to generalize somewhat, band-limited 
noise, limited to the test systems design data-handling capabilities, seems 
to be a reasonable check — primarily because it ’’exercises” the entire design 
spectrum simultaneously. 

For the sake of discussion, one may first consider the two basic types 
of data with which to load a system to. be tested, narrowband and wideband. 

By narrowband, one typically uses a sine wave signal (or perhaps sum of sine 
waves) and measures the delay and distortion of the system output compared 
to the input to get the induced distortion. Wideband data is meant to be data 
which will simultaneously exercise the system under test throughout its 
entire designed spectral range. It may be generated by a filtered noise source 
to produce the desired spectrum and a near-Gaussian amplitude probability 
density distribution, a sum of sine waves, or a filtered PN sequence. 

Waveform distortion in a telemetry system or subsystem may be 
defined as the minimum mean-squared difference between the input signal 
and the output signal. Constant gain or delay between the input and output is 
not considered as error; therefore, an efficient means of eliminating these 
factors from the error calculation must be devised. Mathematically, the 
error is defined as 


cr 2 (t) = ^ / [e.(t - r) - Ge 0 (t) - Gn 0 (t )] 2 dt 
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where e. (t - r) is a delayed replica of the input signal, e 0 (t) is the output 

signal, and n Q (t) is the output noise. To determine the waveform distortion, 
the delay and the gain G must be that of the system under test. 

Two methods, and the basic reason for this appendix, for measuring 
waveform distortion are presented. The first approach follows the general 
ideas of Shepertycki and uses analog techniques for computing the error. The 
second approach employs digital techniques for determining the error. 


Analog Method 

Figure A-l outlines a method for measuring signal distortion first 
proposed by Shepertycki [7] . The test signal is generated as a discrete 
pseudorandom digital waveform clocked through a series of shift register 
elements so that each element has an output delayed by exactly one clock 
period relative to the preceding element. The test signal, simulating the 
wideband data, is formed by passing the binary output from stage one of the 
pseudorandom signal source through a selected lowpass filter, called the data 
filter. A delayed replica of the pseudorandom signal is applied to a second 
data filter, identical to the first, and then compared with the communication 
system output in a differential amplifier. The two pick-off points for generating 
e. (t) and e. (t - r) are placed such that the delay between them is 



Figure A-l. Analog technique for measuring signal distortion. 
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approximately equal to the delay through the communication system. The 
clock frequency is then varied for fine adjustments of the delay. 

Both delay and gain are adjusted to minimize the differential amplifier 
output. This minimal output represents the waveform distortion introduced 
by the communication system and any residual errors contributed by the test 
apparatus. The degree of residual errors is easily evaluated by repeating 
the measurement with the communication system removed. 


Application of the analog technique requires careful selection of clock 
frequency and the pseudorandom signal generator. In general, the clock 
frequency should be made at least 10 times greater than the 3 dB bandwidth 
of the data filter. This insures that the data spectrum shape is a function 
primarily of the frequency characteristics of the data filter. The number of 
stages in the PN sequence generator and delay circuitry should be chosen so 
that the total time delay available in the generator is at least equal to the 
delay in the system under, test. 

The equipment involved with the analog technique has the advantage 
of being compact and relatively inexpensive. On the other hand, some 
undesirable characteristics associated with this technique are (l) the data 
filters must be matched very accurately, which in practice has been difficult 
to achieve, and (2) the analog computing devices are inherently nosiy and 
subject to drift and frequency distortions. 

Some additional problems may be encountered in maintaining the 
desired signal characteristics. Provided the bandwidth of the analog shaping 
filter is much less than the clock frequency, the test signal possesses a flat 
frequency spectrum out to the bandwidth of the filter and (supposedly) has a 
Gaussian amplitude probability density distribution. Unfortunately, in many 
applications involving this method of noise generation, the output becomes 
skewed for large clock frequency/data filter bandwidth ratios. 

Digital Method 

With some complexity, sample-and-hold circuits can be used to pro- 
vide the necessary time delay for measuring signal distortion. A block 
diagram of a unit employing these so-called digital techniques is shown in 
Figure A -2. In this system, time delay is adjusted continuously by variation 
of a clock frequency. The clock ,is used to drive an input scanner (time- 
division multiplexer) and an analog-to-digital (A/D) converter. The input 
scanner is composed of ring counter elements that drive analog gates leading 




Figure A -2. Digital technique for measuring signal distortion. 


to the sample -and -hold circuitry of a digital-to-analog (D/A) converter. The 
time delay between any two ring-counter elements is an integral number of 
clock periods so that a direct relationship between time delay and clock period 
is established. Thus, the clock period is easily converted to units of time 
delay. 

In the unit outlined in Figure A-2, (to) is a sampled version of 

the communication system output and E (t 0 - r) is the corresponding delayed 

version of the input signal. Both outputs appear in the form of a 10-bit PCM 
serial wavetrain. The Telemetry Data Analyzer 1 synchronizes to this wave- 
train, separates the two data samples, and performs the necessary operations 
to determine the error. The number of samples taken for a calculation 
depends on the required accuracy and may be selected as desired. 

Comparison of Results 

It should be noted that any analog device may be inserted into the 
"system” section of the test configurations in Figures A-l and A-2 and its 
error determined. To establish some correlation between these suggested 


1. The Telemetry Data Analyzer is special purpose equipment designed for 
Astrionics Laboratory to provide, among other things, a means of rapidly 
performing sampled rms end-to-end error evaluations. 
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error measurement techniques and analytical results, the error resulting 
from a 6 -pole Butterworth lowpass filter was measured. With the analog 
technique, the pseudorandom output from a 20-stage generator was passed 
through a 3 -pole Butterworth lowpass filter to form the data source. With 
the digital technique, white noise shaped with a 3-pole Butterworth lowpass 
filter provided the data source. In Figure A-3, rms error is expressed as 
a percent of rms input signal strength. It is seen that theoretical and experi- 
mental results agree closely, giving a great degree of confidence in both error 
measurement techniques. 

The techniques described have been successfully used to measure 
single-channel and multi-channel distortion CBW/FM telemetry systems. 

As expected, the results offered a high degree of reproducibility, leading 
to substantially the same results with each successive measurement. 




Figure A -3. Distortion of third-order band-limited data by a 
six -pole Butte rworth lowpass filter. 
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APPENDIX B 

EVALUATION OF A MATCHED LOWPASS FILTER SET 

This appendix gives the results of an evaluation of a pair of matched 
lowpass filters (manufactured by Kappa Networks, Carteret, New Jersey). 
The purpose of the evaluation was to determine if the quality of the filter 
match was sufficient to qualify the filters for application in the pseudonoise 
test set. 

Figure B-l gives a standard schematic for lowpass filters of the type 
evaluated. The filter specifications were as follows: 

Filter Order — Third 

Type Response — Lowpass, Butte rworth 

Load Impedance — Rj_ = 5000 ohms 

Match Quality — 0. 05 percent 



Figure B-2 indicates the procedure used by Kappa in checking the 
match quality of the filters during tuning. With a sinusoidal input signal, 
the filter set was checked throughout its frequency range of operation for a 
normalized difference signal or 0.0005. The frequency range of operation 
was taken to be about three times the dli cutoff frequency. The assumption 
was made that if 



fj (co) - £2 (co) <■ 0.0005 f t (co) 


(B-l) 


for all, where (co) is filter number 1 output and ^ (co) is filter number 
2 output, then 


f [ f 4 (co) - f 2 (co)] dco ^ 0.0005 f fj (co) dco 
0 0 


(B-2) 


or 


f [f 1 (co) - f 2 (co) ] dco 

_0 

CO 

f fj (co) dco 
0 


< 0 . 005 , 


(B-3) 


and the filter set would meet the desired match quality specification. 



BATTERY 
OPERATED 
AS AN AC 
DIFF. 

VOLTMETER 
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Figure B-2. Kappa's test procedure. 





To check the final match quality of the delivered filter set, an input 
drive signal with an approximately uniform spectrum was used. 

Figure B-3 gives the schematic of the test circuit. The drive signal 
was a pseudorandom digital sequence with a bit rate of 6 kHz. The sequence 
length was 1023 bits, and this resulted in an input data spectrum with a 
harmonic component spacing of 5.86 Hz. 

o 

The filter test procedure was as follows: 

1. Close switch Sj and adjust the gain control to confirm the ability 
of the difference amplifier to null to better than 0.01 percent normalized. 

2. Open switch Sj and adjust the gain to the best null. Measure the 
rms value of the difference signal V^. 

3. Close switch s 2 and measure the rms value of the reference signal 

V . 
r 

4. Calculate match quality in percent, (V^/V^ 100 percent. 

The measured match quality was not quite as good as the 0. 05 percent 
specified and was found to vary slightly with the magnitude of the input signal. 
Typical measurements were: 


v V , 

r d 

Match Quality 

0.95 0.001 

0.11% 

0.66 0.0008 

0.12% 


Even though the measured match quality was not as good as the speci- 
fied value of 0. 05 percent, it may still be within the limits allowable in the 
pseudonoise test set. 
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FILTER 



Figure B-3. Filter test scheme 
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